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ABSTRACT 

Duration and period of transits in extrasolar planetary systems can exhibit long-term variations for a 
variety of reasons. Here we investigate how systemic proper motion, which steadily re-orients planetary 
orbit with respect to our line of sight, affects the timing of transits. We find that in a typical system 
with a period of several days proper motion at the level of 100 mas yr _1 makes transit duration vary at 
a rate ~ 10 — 100 ms yr _1 . In some isolated systems this variation is at the measurable level (can be as 
high as 0.6 s yr" 1 for GJ436) and may exceed all other transit timing contributions (due to the general 
relativity, stellar quadrupole, etc.). In addition, proper motion causes evolution of the observed orbital 
period via the Shklovskii effect at a rate > 10 /is yr^ 1 for the nearby transiting systems (0.26 ms 
yr^ 1 in GJ436), which in some cases exceeds all other contributions to P. Earth's motion around 
the Sun gives rise to additional periodic timing signal (even for systems with zero intrinsic proper 
motion) allowing a full determination of the spatial orientation of the planetary orbit. Unlike most 
other timing effects the proper motion signatures persist even in systems with zero eccentricity and 
get stronger as the planetary period increases. They should be the dominant cause of transit timing 
variations in isolated wide separation (periods of months) systems that will be sought by Kepler. 

Subject headings: astrometry — celestial mechanics — planetary systems — eclipses 



1. INTRODUCTION. 

Planetary transits have provided us with a unique op- 
portunity to get a handle on the physical properties of 
the extrasolar planets such as their radii and densities. 
Recently it has been suggested (Miralda-Escude 2002; 
Heyl & Gladman 2007) that precision timing of the mo- 
ments at which transits occur can give us additional in- 
formation about the transiting systems. Various physical 
effects cause orbit of the planet precess in space leading 
to the changes in transit geometry, which can be mea- 
sured through the timing of transits. Among these effects 
are the general relativistic precession of the orbit, gravi- 
tational influence of other planets in the system or com- 
panion stars, torques due to the spin-induced quadrupole 
moment of the star and due to the tidal deformations of 
both the star and the planet (Miralda-Escude 2002; Heyl 
& Gladman 2007; Ribas et al. 2008; Pal & Kocsis 2008; 
Jordan & Bakos 2008). 

Another obvious reason for the re-orientation of the 
planetary orbit with respect to observer at Earth is the 
proper motion of the exoplanetary system with respect 
to the Solar System. Stars in the Solar neighborhood 
move at velocities of tens of km s _1 and some of them 
exhibit proper motion at the level of 1 mas yr _1 . Also, 
the distance to stars constantly changes as a result of 
their relative motion with respect to the Sun and this 
affects transit timing because of the finite speed of light. 
At some level proper motion is a characteristic of any 
star, including those with transiting exoplanets, and it is 
thus important to understand its implications for transit 
timing. 

Proper motion is well known to be important in the 
timing of isolated and binary radio pulsars (Shklovskii 
1970; Kopeikin 1996). In these systems proper motion 
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affects the pulsar spin and orbital periods through the 
so-called Shklovskii effect (Shklovskii 1970) while the re- 
orientation of the binary orbit can be (and has been) 
measured via the variation of the projected size of the 
orbit (Kopeikin 1996; Arzoumanian et al. 1996). Pul- 
sar acceleration in external gravitational field can also 
be important especially for pulsars in globular clusters 
(Edwards et al. 2006). 

Of course, there are significant differences between the 
timing of pulsars and of planetary transits: accuracy 
with which some millisecond pulsars can be timed is 
at the /xs level (Manchester 2008) while a single plane- 
tary transit can only be timed to several seconds at best 
(Knutson et al. 2007). Also, the whole idea of timing 
is different in the two cases: for binary pulsars one is 
usually able to trace the whole orbit of the neutron star 
in time domain while in the case of planetary transits 
only two narrow time windows — primary and secondary 
transits — are available to play with. Nevertheless, some 
of the ideas developed in pulsar timing may be applied 
to the timing of planetary transits. 

Previously, Kopeikin & Ozernoy (1999) have utilized 
a post-newtonian relativistic approach for the precision 
Doppler measurements of the binary star orbits and dis- 
cussed some of the relevant effects of the proper motion. 
Here we aim at investigating the role of the proper mo- 
tion in timing of planetary transit duration and period 
in extrasolar planetary systems. We lay out the basics of 
the orbital element evolution due to the proper motion 
in $2] . We discuss the evolution of the transit duration 
in fj3] and the evolution of the period between transits 
in iJU Comparison with other transit timing effects and 
application to real systems can be found in SJS1 

2. EFFECT OF PROPER MOTION. 

To quantitatively evaluate the effect of stellar proper 
motion on the timing of planetary transits let us consider 
a planet in orbit around a star with period P, semimajor 



axis a, and eccentricity e. We introduce a unit vector n 
pointing from the observer at Earth to the barycenter of 
the transiting system. Vector n varies in time because of 
the linear motion of the binary: 
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(i) 



where fj, is the proper motion in the plane of the sky. Ori- 
entation of the binary in space is fully determined by the 
unit vector 1 parallel to the orbital angular momentum L 
of the binary (i.e. 1 is perpendicular to the orbital plane) 
and the unit vector g pointing from the prime focus of the 
planetary orbit towards its pericenter. We assume L to 
be constant thus neglecting possibility of tidal coupling 
between L and planetary and stellar spins, and gravita- 
tional effects of any companions. We also assume that 
orientation of the orbital ellipse in space is fixed, i.e. g 
is constant too. In doing this we disregard precession 
of the planetary orbit caused by the general relativity, 
stellar oblateness, and so on. We can do this because 
observed changes of the orbital configuration caused by 
different physical mechanisms add up linearly and here 
we want to concentrate on just one of them. 

Orbital plane crosses the plane of the sky along the line 
of nodes and we introduce vector m = (lxn)/ sin i along 
this line (|m| = 1), where i is the observed inclination 
of the planetary orbit given by sini = |n x 1|. If lo is 
the angle between m and g in the direction of planetary 
motion — the argument of pericentre — then at any 
moment of time 

cosw,, , sinw. ,, ,.. . . 

g=— lxn-— n-ln-1. (2) 

sin i sin i 

Differentiating relation cos i — (n • 1) with respect to 
time we find using equation ([I]) 



(m-i) 



-jjl cos /3, 



(3) 



where /i = |/x| and (3 is the angle in the plane of the sky 
between [i and vector 1 — n(l ■ n) — the projection of 1 
on the sky plane. Differentiating with respect to time 
relation cosw = (g • m) = (g • (1 x n))/sini and using 
equations d), ©, and (J3]) we find (Kopeikin 1996) 
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Equations ([3]) and f4| fully determine evolution of the 
observed orientation of planetary orbit in space caused 
by the stellar proper motion. 

3. VARIATION OF THE TRANSIT DURATION. 

Planet transit is characterized by an impact parame- 
ter p — r tr cosi/R* — minimum separation between the 
planetary trajectory and the stellar disk center projected 
onto the plane of the sky, in units of stellar radius R+. 
Here r tr is the value of the spatial separation r between 
the planet and the center of the star at transit midpoint 
- moment of time when the projected separation be- 
tween the planet and the center of the stellar disk is 
minimized. In general 



q(l-e 2 ) 

1 + e cos / ' 



(5) 



where / is the true anomaly counted from the line of 
apsides. Transit midpoint occurs at / = tt/2 — lo, so 
that 

a cosifl — e 2 ) 

P= " -■ 6) 

it* 1 + e sin lo 

Clearly, for the transit to occur one needs p < 1 + R P /R+, 
where R p is the planetary radius, which translates into 

R+ + R„ 1 + e sin lo 

cos«< y -— — . (7 

a (1 — e z ) 

Transit duration calculated as the time between the 
crossings of the edge of the stellar disk by the center 3 of 
the planetary disk is (see e.g. Tingley & Sackett 2005) 



T tr — 2 
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where ify, tr = na(l+esinu)/(l — e 2 ) 1 ' 2 is the value of the 
azimuthal (transverse) component of planetary velocity 
at the transit midpoint and n = 2tt/P is the planetary 
mean motion. In deriving equation (|8|) we have neglected 
the curvature of projected planetary trajectory and the 
variation of planetary speed during the transit — this 
introduces only a small error. 

Given that p and lo in equation J8]) evolve as a result 
of stellar proper motion it is obvious that T tr would not 
remain constant. Differentiating expression (|8|) with re- 
spect to time one finds 

rr Ttr 
J-tr — —~ : 
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where 
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(10) 



R* i - r 

In equation ([9]) the first term in brackets describes the 
variation of T tr caused by the change of v v-tr due to the 
precession of the orbital ellipse while the second and the 
third terms embody the variation of transit geometry 
(change of impact parameter p) caused by the change of 
the inclination of the orbital plane and the precession of 
the orbital ellipse respectively. Third term is normally 
much smaller than the second one because cosi < 1 in 
transiting systems. Note that u> affects T tr only if the 
planetary orbit is eccentric, while i causes variation of 
T tr even for circular orbits. 

Expression for T tr caused by the proper motion can be 
written with the aid of equations ©, ([!]), and ^ as 

T tr p sin (3 
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(11) 



sinj " \tan/3 sini 1 + esino;. 

This equation explicitly shows how T tr varies as a func- 
tion of the absolute value of the stellar proper motion p 
and the orientation of p with respect to the projection 
of the orbital angular momentum onto the plane of the 
sky — angle (3. 

3 This definition differs from that usually adopted in the litera- 
ture (which assumes that transit lasts while stellar and planetary 
disks have at least some overlap) but this does not affect our results 
significantly. 



4. VARIATION OF THE ORBITAL PERIOD. 

Precession of planetary orbit makes observed orbital 
period of the planet P b s different from the true orbital 
period P. Indeed, suppose that we try to determine the 
period of transiting planet by measuring the time be- 
tween the successive inferior conjunctions of the planet. 
Every orbital period precession at a uniform rate u turns 
the orbit by an angle A<p = Plo with respect to our 
line of sight which gets reflected in the length of the 
time interval between successive conjunctions. The ex- 
tra time it takes a planet to cover this additional angle 
is AP^ = Aip/ip, where tp = v v j r /r tr is the angular fre- 
quency of the planet at the point of conjunction. Using 
equation ((5|) and expression for v Vl tr we can write the 
deviation of observed planetary period from the true one 
as 



AP, 



(1-e 2 ) 
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(12) 



Note that AP U is nonzero even in the case of circular 
orbits in agreement 4 with Kopeikin (1996) and Heyl & 
Gladman (2007). In general Co in equation (TT2"]) is given 
by 



LO = Wu +LUGR +&S, 



(13) 



where different terms on the right hand side represent ap- 
parent orbital precession caused by the systemic proper 
motion, general relativity and the quadrupole moment of 
the central star correspondingly (other sources of orbital 
precession, e.g. tidal deformations of star and planet 
(Jordan & Bakos 2008) have been neglected here for sim- 
plicity). As a result, AP U = AP U ^ + AP UjGR + AP UjS , 
where AP U u, AP U gr and Ai^^ are found by substi- 



tuting Up, ujgr, and Cos for Co in equation (|12|). 

Apart from AP U which owes its existence to the appar- 
ent re- orientation of the planetary orbit there is another 
contribution to P bs related to the systemic motion: the 
distance to the planetary system changes, which because 
of the finite speed of light gives rise to a special relativis- 
tic contribution AP re i given by 



AP, 



rel 



p 



c 



(14) 



where v r is the line-of-sight velocity of the system (posi- 
tive for systems moving away from us). Thus, in general 

P obs = P + AP^ + AP rel . (15) 

One cannot, of course, measure AP directly since the 
true orbital period of the system is not known a priori. 
However, one might try to measure the variation of P 
over an extended period of time. With equations (fT2|) . 
(TT4"|) , and (fT5|) one can easily show that 



P = P, 



P 
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where 
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(16) 
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4 By contrast expressions for AP U derived in Miralda-Escude 
(2002) and Pal & Kocsis (2008) vanish in the limit e -> 0. 



with D being the distance to the planetary system and 
Vt = fiD being its transverse velocity. The timing con- 
tribution Pshk, which to the best of our knowledge has 
never been highlighted in the context of planetary tran- 
sit timing, is identical to the so-called Shklovskii effect 
well known from pulsar timing (Shklovskii 1970): radial 
motion of the system changes the observed orbital pe- 
riod via the Doppler effect but the radial component of 
the velocity (and the Doppler factor) varies if there is a 
non zero transverse component of the systemic velocity, 
leading to non zero P. This contribution to P is always 
positive since spatial motion of the planetary system al- 
ways increases v r . In $5] we demonstrate that in many 
transiting systems Shklovskii effect dominates P. 

Equation ([17]) generalizes expressions for P^ obtained 
by Heyl & Gladman (2007) and Pal & Kocsis (2008) 
by including the term proportional to Co. One expects 
Co <C {to) 2 for precession caused by the general relativ- 
ity and the stellar quadrupole in which case equation 
(fT?) reduces to the expression derived by other authors. 
However, in the case of apparent precession caused by 
the systemic motion one can easily show using results of 



that 



fj, 2 ~ (tip) , so that all terms in equation 



(fT7|) for P W)At must be retained. In general, P w <C T tr 
because T tr is a linear function of the small parameter 
Plu while P u is quadratic. 

5. DISCUSSION. 

Here we compare the effects caused by the proper mo- 
tion with other timing contributions and discuss their 
observability in different types of systems. 

As a fiducial system we will take a star located 100 pc 
away from the Sun and moving with transverse velocity 
30 km s _1 . Such a system has proper motion /i « 60 mas 
yr _1 resulting in « M , w M rs 2 x 10~ 7 yr _1 for = 45° and 
i = 90°. Timescale on which planetary orbit changes its 
orientation is ~ /i -1 ~5x 10 6 yrs. We can compare lu^ 
to the general relativistic periastron precession rate 
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and to the rate of precession due to the rotation-induced 
stellar quadrupole (Miralda-Escude 2002) 

7/2 
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where Ji is the dimensionless measure of the stellar 
quadrupole moment (its typical value for the Solar type 
stars is J2 ~ 10~ 6 ) and we took M* = M and R+ = i? Q . 
These estimates clearly indicate that for Solar type 
stars with short period (P = 3 — 4 days) planets logr. ^> 
^s 3> ifi, LOft. Plugging expression (fi9|) into equation (fT7| 
we find that a planetary system with M* = M®, e = 0.1 
and lo = 45° should exhibit 



P 



367recosw 



u,GR' 
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(21) 



us yr 



Given that u^ <§C Cjgr it is clear that P W , M <C Pu,gr so 
that the re-orientation of the planetary orbit caused by 
the stellar proper motion does not noticeably affect P u 
(the same is true for the precession caused by the stellar 
quadrupole since u>s <C ^>gr)- 

However, this does not mean that one can just ignore 
the effect of the proper motion on P: proper motion also 
affects P via the Shklovskii effect and the magnitude of 
this contribution 
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may be comparable to P u ,gr- Clearly, Pshk can be quite 
important even for tight, eccentric systems for which one 
would normally expect P^.gr to dominate. 

One also has to keep in mind that the majority of short 
period transiting systems have eccentricities consistent 
with zero. In such systems with circular orbits P^^gr 
and Puj.s vanish (remember that u>gr,&S ~ 0) leaving 
Shklovskii effect as the only source of non zero P w at the 
level of tens of /us per year. In Table [TJ we have summa- 
rized the properties of observed transiting systems (sup- 
plemented with two artificial systems Sys-1 and Sys-2 
with the goal of illustrating transit timing effects in long 
period systems) in which proper motion effects are par- 
ticularly pronounced (namely, max|T t ,. iA ,| > 10 ms yr^ 1 ), 
while in Table [2] we display the values of various tim- 
ing contributions in these systems, including P^^gr an d 
Pshk- From Table [2] one can see that in some nearby 
high-proper motion systems like GJ436 Pshk is a good 
fraction of ms yr -1 . Such high rate of period change sig- 
nificantly exceeds P u ,gr and may in principle be measur- 
able on a time scale of tens of years assuming observing 
parameters typical for the Kepler photometric mission 
(Miralda-Escude 2002; Jordan & Bakos 2008). 

Variation of the transit duration T tr presents another 
way of detecting proper motion effects in isolated star- 
planet systems, as described in |J31 Assuming that all an- 
gle dependent factors in equation (fTTj) are of order unity 
one finds 
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where in evaluating g we have assumed p = 0.5-R*. Thus, 
a typical nearby exoplanetary system indeed exhibits 
Ttr. ii 3> P ■ A specific value of T tr ^ for a particular 
exoplanetary system depends not only on fi but also (si- 
nusoidally) on the angle (3 between \x and the line of 
nodes. The maximum possible value of T trfi for several 
representative systems can be found in Table [2] 

At the same time, for M* = M Q , e = 0.1 and u> = 45° 
one finds from equation ([9]) the following value of T tr due 
to the general relativity: 



f tr ,GR pa 240 
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This is not much larger than T tr .^ and in some high 
proper motion systems T tr .^ may even dominate. The 
best example is GJ436: because system is very compact 
general relativity provides T tr ,GR ~ —0.2 s yr" 1 but the 
very high proper motion of the system (/x ~ 1.2 mas 
yr" 1 ) gives rise to max|Ti rj/i | w 0.6 s yr . Thus, in 
general one cannot simply ascribe all T tr measured in ec- 
centric systems to the general relativity — some fraction 
of T tr can also be contributed by the proper motion. In 
systems with circular orbits T tr .GR = 0. 

Given that lu^ <C ojgr the magnitude of the effect of 
the proper motion on T tr may seem disproportionately 
large compared to T tr ,GR- The reason for this lies in the 
amplifying factor g in equation §§§ which propagates into 
T{r,/i> see equation (fTTj) . According to the equation © 
the magnitude of g is determined by the transit impact 
parameter p and the ratio a/ R± which is usually quite 
large, ~ 10 even for rather short period (P = 3 — 4 d) 
systems. For grazing transits (such as those occurring in 
GJ436, see Table[2]), when 1 — p -C 1, g gets additionally 
boosted up 5 because then T tr becomes a very sensitive 
function of p and i, see Ribas et al. (2008). At the same 
time factor g ^> 1 does not greatly affect T tr ,GR since for 
precession induced by the general relativity i = and 
g enters the expression for T tr ,GR only in combination 
g cos i while cosi <C 1 in transiting systems (cosi < i?*/a 
so that g cosi ~ 1, see equations ^ and (fTTJj) ) . This 

explains why T tTt GR ~ T tr ^ even though Cjgr ~^> w^. 

Inclination of the planetary orbit with respect to our 
line of sight may change also because of the spin in- 
duced quadrupole if the stellar spin axis is misaligned 
with the orbital angular momentum vector. The spin in- 
duced T tr ,s is amplified by factor g in a way analogous 
to the amplification of T tr .^- Given that in some systems 
u>s can be 1-2 orders of magnitude larger than io^ (sec 
Table [2] where d>s is computed for J-2 = 10~ 6 ) one may 
expect Ttr,s 3> T tr ^ in these systems. However, in re- 
ality it will often be the case that is ^C ws since it can 
be demonstrated that is = Cws sin A (Lai et al. 1995), 
where C ~ 1 is the angle-dependent factor and A is the 
angle between the stellar spin axis and the orbital an- 
gular momentum vector. Misalignment angle A has been 
measured in several systems via the Rossiter-McLaughlm 
effect (Rossiter 1924; McLaughlin 1924) and in the ma- 
jority of measured cases A is close to zero, as expected 
from the planet formation theories. Among the systems 
in Table [TJ for which A has been measured this angle was 
found to be small in HD189733 (A = 1.4° ± 1.1°, Winn 
et al. 2006) and HAT-P-1 (A = 3.7° ± 2.1°, Johnson et 
al. 2008) while in HD17156 misalignment may be sig- 
nificant (A = 62° ± 25°, Narita et al. 2008), although 
Cochran et al. (2008) have found A = 9.4° ± 9.3° in this 
system. In HD189733 is end up being <C i M so that T tr ,s 
likely makes negligible contribution to T tr which should 
be dominated by the proper motion. In HAT-P-1 we find 

is ~ V anci T tr,s ~ max|T tr ^| with f tr ,GR providing a 
non-negligible contribution. Finally, in HD17156, if we 

5 Thus, monitoring systems with large p tends to increase the 
chances of measuring Tt r . 



adopt a larger value of A found by Narita et al. (2008), 
sin A ~ 1 but is is still comparable to i^ because the 
semimajor axis of the system is quite large which greatly 
reduces lus- As a result, Ttr.s ~ max|Tt ri/i | in this sys- 
tem and both are somewhat smaller than T tr .GR- Thus, 
at least in the systems presented in Table [T] the spin- 
induced quadrupole orbital precession does not strongly 
exceed the proper motion effects in timing of transit du- 
ration. 

Note that the tidal deformations induced on the star 
and the planet by each other affect T tr in a way different 
from that of the spin induced quadrupole — similar to 
the general relativity the tidal bulges do not generate 
non zero i. Given that tidal u> is typically smaller than 
ugr (Jordan & Bakos 2008) we may conclude that tidally 
induced T tr is lower than T tr .GR and is thus < Ttr.fi- 

It is obvious from the preceding discussion that the 
proper motion can have an appreciable (if not dominant 
in some cases) effect on transit timing in the short pe- 
riod systems. This statement becomes much more robust 
when we go to systems with wider separations. It is obvi- 
ous from equations 1(19]). (j20|) . (|2Tj). and J24} that oj G r, 
ibs, and all contributions to P and T tr caused by the 
effects of the general relativity, stellar quadrupole and 
tidal deformations are rapidly decreasing functions of a. 
At the same time, tl> M , i^ are independent of a while both 
Pshk and T tr increase quite rapidly with a, see equa- 
tions (|22|) and (f23|) . This means that the proper motion 
should completely dominate transit timing variations in 
isolated (i.e. containing no other planets) wide separa- 
tion systems. For example, a transiting planet in a 30 d 
orbit around a Solar type star would exhibit Pu.gr ~ 0.6 
/is yr -1 and T tr .GR ~ 7.5 ms yr -1 if e = 0.1 and T tr = 4 
hr. If this system is located 100 pc away from the Sun 
and has proper motion /i = 100 mas yr" 1 then one finds 
Pshk ~ 200 /js yr -1 and T tr ^ « 200 ms yr -1 , so that 

both Puj.gr < Pshk and T tr .GR < T tr ^. To additionally 
illustrate the importance of proper motion for wide sepa- 
ration systems we introduce two artificial systems (Sys-1 
and Sys-2) in Table[T]and calculate their timing parame- 
ters in Table HJ Such wide separation systems are one of 
the primary goals of photometric missions like Kepler. It 
is clear that if such systems are found to exhibit transit 
timing variations then these variations must be caused 
by the systemic proper motion, provided that the influ- 
ence of possible additional companions is proven to be 
negligible. In this case according to equation (jTTJ) the 
measurement of T tr would serve as a measurement of an- 
gle /? giving us information on the full three dimensional 
orientation of the transiting system. However, one must 
remember that if additional planets in external orbits are 
present in these systems then their influence may not be 
disregarded (Ribas et al. 2008) since their effect on tran- 
sit timing grows with a faster than T tr ^ and Pshk do. 

Note that the rapid increase of Pshk and T tr .^ with a 
does not immediately imply that their actual detection is 
facilitated as P increases. Even though the timing signal 
increases with P, the number of transits, which deter- 
mines the timing error, decreases as P~ l for a given time 
interval over which the system is being monitored. Us- 
ing the results of Ford et al. (2008) and Heyl & Gladman 



(2007) on transit timing precision we find that the time 
lip one needs to monitor the transiting system to detect 
P induced by the proper motion scales as 6 lip ex P 4 / 15 ; 
i.e. it increases with P but not very rapidly: it takes 3.6 
times longer for P = 1 yr system to get the same S/N for 
Pshk due to proper motion as for the 3 d system. The 
uncertainty in T tr decreases with the number of observed 
transits slower than the uncertainty in P. As a result, 
the time II tr one needs to monitor the transiting system 
to detect T tr caused by the proper motion decreases with 
P as 7 Utr oc p -2 / 9 . Thus, it is easier to measure T tr . M 
in wide separation transiting systems. 

Based on the results of Ford et al. (2008) it was esti- 
mated by Jordan & Bakos (2008) that the Kepler mission 
should be able to achieve a timing precision of ~ 1.5 s 
in 1 yr observation of a 12th magnitude Solar type star 
transited by a Jupiter-like planet with P = 5 d. One can 
easily deduce from this that a 3a detection of T tr ^ = 100 
ms yr -1 (which is not unreasonable for proper motion) 
should take H tr « 10 yr of observations. Measurement 
error of P drops very rapidly with time but it would still 
take lip w 70 yr to achieve a 3a detection of P = 100 
fis yr -1 caused by the Shklovskii effect. Thus, while one 
might hope to measure Ttr,p, in some nearby, high-proper 
motion systems (like GJ436) on timescale of ~ 10 yr, the 
measurement of P would likely require next generation 
facilities with photometric precision much higher than 
that of the Kepler mission. 

In systems with low proper motion T tr ^ can be viewed 
as an irreducible systematic uncertainty to which quanti- 
ties like Ttr.GR can be measured. This is because even if 
/i is known precisely one still does not know a priori the 
angle /3 (but see below) which determines T tr ^. Thus, 
proper motion limits to some extent our ability to in- 
terpret the measurement of T tr in terms of the physical 
parameters of the system (e.g. J2, etc.). Measurement of 
P does not suffer from this uncertainty since Pshk is in- 
dependent of and can thus be fully accounted for once 
\i is known from astrometric measurements. 

Previous discussion has implicitly assumed that a tran- 
siting system moves at a constant speed with respect to 
observer. In reality observer is located at Earth, which 
orbits the Sun. This to some extent complicates the anal- 
ysis of the transit timing data since the relative velocity 
between the observed system and the Earth is a function 
of time with a 1 yr period. But given that the Earth-Sun 
motion is known its effect on transit timing can be easily 
accounted for: terrestrial orbital motion generates pe- 
riodic apparent proper motion of any transiting system 
(even if it has no intrinsic proper motion) with respect 
to observer at Earth. The maximum amplitude of this 
apparent proper motion is /ip ~ ve/D w 60 mas yr" 1 , 
where tp w 30 km s" 1 . According to equations (JTTJ) 
and (|23|) \xe gives rise to periodically varying T tr ^ with 
an amplitude dependent on the orientation of the orbital 

6 This dependence is found by equating Pshk^p to opi/iV where 
N = Tlp/P is the number of observed transits and op is the un- 
certainty in measurement of Ttr, which is given by eq. (6) of Ford 
et al. (2008). 

7 This follows from equating Ttr.^Tltr to op given by eq. (6) of 
Ford et al. (2008). 



Table 1. Parameters of systems with transiting planets 



System 


/la 
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P 


a 
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U) 


i 


R* 


M* 




(mas yr~ 1 ) 


(mas yr -1 ) 


(pc) 


(d) 


(AU) 




(°) 


(°) 


(R&) 


(Mq) 


HD149026 


-77.12 


53.34 


74.4 ±7.2 


2.87588 


0.0432 





_ 


85.3 


1.368 


1.294 


HD189733 


-2.49 


-250.81 


19.7 ±1.0 


2.21857 


0.0312 





- 


85.58 


0.756 


0.806 


HD17156 


91 


-32.45 


78.24 


21.21725 


0.1594 


0.6717 


121.23 


88.23 


1.47 


1.2 


GJ436 


896.34 


-813.7 


10.2 ±0.2 


2.64385 


0.02872 


0.15 


351 


86.5 


0.464 


0.452 


TrES-4 


-8.1 


-33.0 


485 ± 31 


3.55394 


0.0488 





- 


82.81 


1.816 


1.394 


XO-5 


-32.4 


-24.4 


270 ± 25 


4.187732 


0.0508 





- 


86.8 


1.11 


1.0 


HAT-P-1 


29.3 


-51.0 


155 ± 15 


4.46529 


0.0551 


0.09 


80.7 


86.11 


1.135 


1.133 


WASP-2 


3.0 


-53.1 


157 ±4 


2.152226 


0.0307 





- 


84.81 


0.84 


0.89 


Sys-1 


7.1 


7.1 


300 


100 


0.4218 


0.5 


do 


89.5 


1.0 


1.0 


Sys-2 


7.1 


7.1 


300 


365 


1.0 


0.5 


45 


89.8 


1.0 


1.0 



plane of the transiting system with respect to the eclip- 
tic, potentially providing a method of measuring angle 
/?. The maximum possible value of such timing signal is 
about 



T, 



tr.E 



30 



/< 



a/R^ T t1 



100 mas yr^ 1 10 4 hr 



ms yr 



(25) 



It is also obvious that terrestrial orbital motion produces 
a periodic contribution to P, even if the transiting system 
has zero intrinsic proper motion. 

Such annual variations in T tr and P can arise only as 
a result of the proper motion effects. One can hope to 
measure them by properly combining the data on tran- 
sit duration measured at different orbital phases of the 
Earth. If such variations are detected then this periodic 
part of the timing signal can be used to constrain angle 
f3 allowing one to remove the aforementioned systematic 
uncertainty in measuring other timing contributions. 

6. CONCLUSIONS. 

We investigated the effect of the proper motion on tim- 
ing of transiting exoplanetary systems in which the grav- 
itational effect of other possible companions can be ne- 
glected. Proper motion re-orients planetary orbit with 



respect to our line of sight and changes the distance 
to the system. Short period transiting systems having 
proper motion at the level of 100 mas yr^ 1 should ex- 
hibit variation of the transit duration at the level of 
~ 100 ms yr -1 , which may be comparable to or exceed 
the timing signatures produced by the general relativity 
or stellar quadrupole and which should not be hard to 
detect. Proper motion also causes variation of the ob- 
served orbital period through the Shklovskii effect which 
dominates P for high proper motion systems. Orbital 
motion of the Earth around the Sun gives rise to pe- 
riodically varying transit timing signal even in systems 
having zero intrinsic proper motion. Timing effects in- 
duced by the proper motion become especially important 
in systems with zero eccentricity and in wide separation 
systems with periods longer than a month which should 
be discovered by Kepler. 
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max uj^j. 
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(10- 7 yr- 1 ) (: 


10- 4 yr- 1 ) (: 


LO- 7 yr- 1 ) 


(ms 


yr" 1 ) 


(ms yr— 1 ) 


(MS 


yr" 1 ) 


(m« 
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HD149026 


5.5 


4.5 


7.1 


114.8 




24.0 







12.46 







HD189733 


11.4 


12.2 


7.8 


23.4 




69.3 







18.2 







HD17156 


6.2 


4.7 


0.44 


1.8 




21.2 


-97.2 




102.6 




-2.4 


GJ436 


31.9 


58.7 


4.1 


3.6 




616.6 


-201.1 




261.8 




13.85 


TrES-4 


8.7 


1.6 


5.8 


297.6 




16.7 







13.2 







XO-5 


7.7 


2.0 


3.2 


26.1 




14.9 







12.3 







HAT-P-1 


11.5 


2.9 


3.2 


20.6 




24.8 


9.8 




15.9 




1.7 


WASP-2 


11.3 


2.6 


9.5 


39.9 




15.4 







6.3 







Sys-1 


49.2 


0.48 


0.02 


0.04 




30.7 


6.7 




19.9 




0.3 


Sys-2 


108 


0.48 


0.002 


0.002 




105.1 


0.7 




72.5 




0.05 



